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^ ■ 1 Introduction 

Robustness is considered as an elementary property in the analysis and design of control systems. 
A control system is said to be robust if it retains an assigned degree of stability or performance under 
perturbations. Robust control under structured perturbations has been an active research area in 
recent yearst H 1. 

Robust stability problem for control systems with interval characteristic polynomials was initiated 
by Kharitonov(1978)^, who showed that the Hurwitz stability of a real interval polynomial family 
can be guaranteed by the Hurwitz stability of four prescribed critical vertex polynomials in this family. 
Bartlett, Hollot and Huang extended the result and presented the Edge Theorem^, which says that 
the stability of a polytope of polynomials can be guaranteed by the stability of its one-dimensional 
exposed edge polynomials. Chapellat et al.(1989)^ refined the Kharitonov's Theorem, and proved 
that the robust stability of the unity feedback control system family with an interval plant and a fixed 
controller in forward path can be guaranteed by the stability of 32 prescribed one-dimensional unity 
feedback control systems in this family. 

When considering multivariable control systems, the robust stability problem becomes more com- 
plicated and challenging [^M^l. For Hurwitz stability of interval matrices, Bialas 'proved' that in 
order to guarantee robust stability, it suffices to check all vertex matrices^. Later, it was shown by 
Barmish that Bialas's result was incorrect [5]. Kokame and Mori established a Kharitonov-like result 
on robust Hurwitz stability of interval polynomial matrices^ . Their result is refined by Wang, Wang 
and Yu^. Kamal and Dahleh established some robust stability criteria for MIMO control systems 
with fixed controllers and uncertain plant [ 8 l 



In this paper, we focus on a class of MIMO feedback control systems with their transfer function 
matrices described by 

M{s)=B(s)A{s)+V{s)C(s) (1) 
**(*) = { + + - + bMe$\*$] } (2) 

w = { 4? + dg'- + • • • + 4°^ <f e Kg ) ,aL- ) ] } ( 3 ) 

B(s) = {(M*))„xn> (4) 

2?( a ) = {(d ij ( S )) nxn , € 2>i»} (5) 

where n is a given positive integer, and b^\b\j\d^\d\j^ are given real numbers. And A(s) = 
(o>ij(s)) nxn ,C(s) = ( c ij( s )) nxn are two fixed polynomial matrices, B(s),V(s) are interval polynomial 
matrices and Bij(s),T>ij(s) are interval polynomials, where A(s) = ( a ij( s )) nxn means that the matrix 
is an n x n one and its ij-th entry is Ojj(s). Rewrite Ai(s) as 

M(s) = M n s n + Mn-is"" 1 + • • • + M 

where Mj is an n x n scalar matrix. 
Assumption A: det(M n ) 7^ over .M(s). 

In the sequel, we study the multivariable interval control system whose characteristic polynomial 
matrix M(s) satisfies Assumption A. By exploiting the uncertainty structures, we are able to reduce 
the computational burden in checking robust stability. 

An outline of the paper is as follows. Some useful definitions and lemmas are presented in Section 
2. Section 3 contains our main result on minimal dimension robust stability criterion. An example is 
worked out in Section 4, and some discussions are made in Section 5. Concluding remarks are given 
in Section 6. 



2 Preliminaries 
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Definition 2.1 Given an interval polynomial set IZ(s) = Y^i!=o r i s% i r i G R > r i\i its Kharitonov 
vertex set is defined as TZ V {s) = {ri(s),r2{s),rs(s),r4(s)}, and its Kharitonov edge set is defined as 
n E (s) = {\ ri (s) + (1 - X) rj (s), G {(1, 2), (2,4), (4,3), (3, 1)}, A G [0, 1]}, where 

n(s) 

7*2 («) 
7*3 («) 

r 4 (s) = r% + s + r^s 2 + r^s 3 + r% s 4 H 

Denote the Kharitonov vertex sets and the Kharitonov edge sets of Bij (s), V, L j (s) as (s) , (s) , Bfj (s) 
and Vfj(s), respectively. 

Definition 2.2 A fix polynomial is said to be (Hurwitz) stable if all its roots lie in the strict left 
half of the complex plane. A polynomial matrix is said to be stable if its determinant is stable. 

Lemma 1 (Edge Theorem^) Let $7 be a polytope of polynomials without degree dropping. Then 
the boundary of R(Q) is contained in the root space of the exposed edges of f2, where R(W) = 
{s such that T(s) = for some polynomial T G W} is called the root space of W. 

Lemma 2 A(s) and C(s) are defined as before. Let B(s) = ipij{s)) nxn ,D(s) = (dij(s)) nxn and 

3fei, k 2 G {1, ■ • • ,n}; l = l,---,n 
B(s)A(s) + D(s)C(s) : b lkl (s) G Bf k As), b u {s) G B%(s), I + h; 



A«(s) = < 



d ik2 (s)evi 2 (s), d u (s)evv(s), l + k 2 



(6) 



3k G {1, • • • , n}; I ^ k, 1=1, ■ • • ,n 
B(s)A(s) + D(s)C(s) : M*) x d lk (s) G (jB£(s) x 2>&( a )) U (s£( a ) x 

M*) x d«( S ) € x X#( a )) 



Then, N^\s) is stable for all z G {1, ■ • • , n} if and only if Ng\s) is stable. 



(7) 



Proof (Sufficiency) Suppose that Ng (s) is stable , our aim is to prove that (s) is stable. For any 
N(s) G ATW(s), there exist /ci,fc 2 G {!,-••, n} such that 



6 ifcl (a) G Bg (s) ; 6« (s) e f# (a) , Z ^ *i , i = 1, 
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,n 



(8) 



N{s) 



( Efe ( bife(s)afei(s) + d lk (s)c kl (s) ) 
Efe ( hk(s)a kl (s) + d ifc (s)c fc i(s) ) 



V Efc ( &nfc(s)ajfci(s) + d nfe (s)c fc i(s) 



Efc ( &ifc(s)a fcn (s) + c?ifc(s)cfe„(s) ) N 

Efc ( M s Hn(s) + d ik {s)c kn (s) ) 

Efc ( Kk{s)a kn (s) + d nfe (s)c fcn (s) ) y 



(9) 



By definition, N(s) is stable if and only if detiV(s) is stable. By Laplace formula, expanding the 
determinant along the i-th row, we have 

detiV(s) = Efc (hk(s)a kl (s) + d ik (s)c kl (s)) N a + • • • + £ fe (M s Wn( s ) + d ik (s)c kn (s)) N in , 
where Nij is the algebraic complementary minor of the ij-th entry of N(s), which is independent of 
bn(s), • • • , bi n (s),dn(s), • • • , di n (s). By a simple algebraic manipulation, we have 

detiV( S ) = E"=i M") (ELi a jk (s)N ik ) + £? =1 (ELi c jk (s)N ik ) 

= b ikl (s)Si(s) + d ik , 2 (s)52(s) + 5 3 (s) 

where <5j(s) is a term which is independent of bi kl (s), di k2 (s). Denote 

' b ikl (s) x d lk2 (s) G ( (Bf ki (s) x Pf fe2 ( s )) U (fl&to x Pf fc2 ( S )) ) 

WeBlW, l^h, l = l,---,n (10) 
I da(a)€X>^(a), I + k 2 , l = l,---,n 
By Lemma 1, the stability of each such polytope is equivalent to the stability of its exposed edges. 
These edges are obtained by letting B(s) and D(s) satisfy (10). Thus, detiV(s) is stable for N(s) 
satisfying (8), (9) if and only if detTV(s) is stable for N(s) satisfying (9), (10). Since N^\s) is stable , 
this shows that N(s) is stable. Hence, sufficiency is proved. 
(Necessary) This is obvious since 



3 Robust Stability Criterion 

In this section, we consider the uncertain system family (1). Before presenting our main result, we 
first define some notations. The notation X xY stands for the Cartesian product, which is the set of 
all ordered pairs (x, y), where x G X, y G Y. Let S n denote the set of all bijections of the set {!,■■■ ,n} 



onto itself, and T n be the set of all functions from the set {1, ■ • • , n} onto itself. Take 

BS (.) = { B (., 6B( .)= w^w : j;^,. e *} 
D , w = { DW ^.):^)^w ; j=^..^} 

= B E {s) x 2? B (s) 

I M^^)e^)x#),jVr,(i) 
The following theorem was given by Kamal and Dahleh (1996). 

Proposition 1® .M(s) is stable for all B(s) G S(s), D(s) G P(s) if and only if M(s) is stable for 
all B(s) x D(s) G S^D^. 

As a generalization of the Generalized Kharitonov Theorem, Proposition 1 addresses the robust sta- 
bility of Ai(s) and reduces it to 4 2n (n!) 2 robust stability problems involving 2n uncertain parameters. 
In this paper, we will show that, by further making use of the uncertainty structure information, the 
original robust stability problem can be reduced to (2n)!4 2n (n!) robust stability problems involving 
n uncertain parameters. By Proposition 1 and Lemma 2, we get the main result of this paper. 

Theorem 1 M(s) is stable for all B(s) G B(s),D(s) G T>(s) if and only if M(s) is stable for all 
B(s) x D(s) G BV. 

Proof: By Proposition 1, M(s) is stable for all B(s) G B(s),D(s) G V(s) if and only if M(s) is 
stable for all B(s) x D(s) G B E V E . For any n(s) G M(s) with B(s) x D(s) G BeT> e , fixing the rows 
from the second to the n-th rows and applying Lemma 2 to the first row, we have 

M(s) is stable for all B(s) G B(s),D(s) G V(s) 

if and only if 

M(s) is stable for all B(s) x D(s) G BV X n B E T> E 
where BT>\ C B(s) x P(s) satisfies (10) when i = 1. Repeating the procedure for the remaining rows 
consecutively, we get the conclusions. 

Remark 1 Theorem 1 shows that, the criterion for the robust stability of interval polynomial 
matrices can be greatly simplified. Proposition 1 of Kamal and Dahleh, reduced the 2n 2 -dimensional 
stability verification problem to 2n-dimensional problem, and our result reduces the 2n 2 -dimensional 
stability verification problem to n-dimensional one. 

Remark 2 Theorem 1 establishes a minimal dimension testing criterion for robust stability of 
multivariable control systems. By making use of the uncertainty information of B(s) and T>(s), i.e., 
not only making use of the uncertainty structure information in B(s) and T){s) separately (this is 
exactly what Kamal and Dehleh did), but also making use of the uncertainty structure information 
collectively, we are able to improve the robust stability results given by Kamal and Dehleh. 

Remark 3 Theorem 1 allows us to improve the main results of Kamal and Dehleh. 



4 Illustrative Example 



Recall the twodink planar manipulator considered by Kamal and Dehleh. 
nomial of the closed loop system is 

T(e) = det (m(^)s 3 + K d s 2 + K p s + K r ) 



The characteristic poly- 
(11) 




M(9 d ) 




where 6 d G [0,7r/2] is the joint angle and 

3 + 2cos(0 d ) l + cos(6> d ) 
1 + cos(0 d ) 1 

k d 22 I 

kp22 

Kl2 
K22 

k dll G [6.07 - 6.07e, 6.07 + 6.07e], k dl2 G [2.22 - 2.22e, 2.22 + 2.22e] 

k d21 G [2.22 - 2.22e, 2.22 + 2.22e], k d22 G [1.62 - 1.62e, 1.62 + 1.62e] 

kpn G [6.12- 6.12e, 6.12 + 6.12e], fc p i 2 G [2.24 - 2.24e, 2.24 + 2.24e] 

k p2 i G [2.24 - 2.24e,2.24 + 2.24e], fc p22 G [1.64 - 1.64e, 1.64 + 1.64e] ' ^ ' 

fern G [5.11 - 5.11e, 5.11 + 5.11e], fe rl2 G [1.87 - 1.87e, 1.87 + 1.87e] 

fc r 2i G [1.87 - 1.87e, 1.87 + 1.87e], fc r22 G [1.37 - 1.37e, 1.37 + 1.37e] 
Since for any fixed 9 d G [0, 7r/2], the inertia matrix 

1 (1 + cos(e d )) 

-l + cos(0 d )) 1 
In what follows, we take the uncertain inertia matrix as 

M(s) = B(s)A(s) 

where A(s) = (\ J ) and 

B(s) G B(s) = {(M*))2x2 : h 3 (s) G 

Bn(s)=/, £ 12 (s) = 6 12 s 3 , 612 €[1,2] 
S 22 (s) = s 3 , B 21 (s) = 6 2 is 3 , 62! G [-1, 0] 

Let V(s) = K d s 3 + K p s 2 + then 

V(s) = {(d ij (s)):d ij (s)eV ij (s)}, 

where 

^V/'(^) — ^^ij(^) — k d ijS 4- kpijS -\- k ri j} 

and k d ij,k p ij and fcrij are given in (12). Thus, the characteristic polynomial of the uncertain system 
is 

T(e) = det (B(s)A(s) + D(s)) ; B(s)eB(s), £>(s) G 2?(a) (13) 
where ^4(s) is a fixed matrix and B(s),V(s) are interval polynomial matrices. 
The Kharitonov polynomial vertex sets associated with Bij(s) and T>ij(s) are as follows, 

BXM = {s 3 } = Bn(s), B\ 2 (s) = {s 3 ,2s 3 }, 

= {s 3 ,0}, BUs) = {s 3 } = B 22 (s), 

( A 1 ) _ K U 2 , x L s 4 K L r/ (2) — K u s 2 -I- K u s -I- K L 1 
ij ~ I / 3) - /r L s 2 4- /T L s 4- K U / 4) - K L s 2 4- K U s 4- K U I 

where the upper bound and the lower bound of K di j are denoted as K^-^K^-. And the Kharitonov 



edge sets associated with B(ij) and V(ij) are as follows. 



B 



(s) 



{s 3 } = B n ( S ), 

{As 3 + (l-A)-2s 3 :Ae [0, 1]}, 
{-As 3 : AG [0,1]}, 
M = B 22 (s) 




A)i 2) , 



Af + 
Xdp + 




By Proposition 1 given by Kamal and Dehleh, for any e, T(e) is robustly stable if and only if T(e) is 
stable for all B(s) x Z?(s) € BeV e , where 

BeVe = U 1 UU 2 



Bu(s) 



Bu 

°2l 



B 22 (s) 



fif 2 ( 



B 



22 



Vf 2 {s) 

vy 2 i 
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» 
*0 



So, the robust stability of T(e) is reduced to 2-4 6 robust stability problems involving four parameters. 
In this case, the testing set is a four-dimensional set. In the sequel, we apply Theorem 1 to simplify 
this problem further. That is to say, T(e) is robustly stable if and only if T(e) is stable for all 
B(s) x D(s) € BV, where 

Vi u v 2 u v 3 u v 4 u v 5 u v 6 u y 7 

Bf 2 («) 
fi 22 (s) 



Vo 
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Clearly, the robust stability of T(e) is reduced to 7 • 4 6 robust stability problems involving two 
parameters. In this case, the minimal testing set is two-dimensional. 

Remark 4: A test with one uncertain parameter could involve infinite tests, e.g., an edge contains 
infinite points. The significance of our result is that it is the minimal dimension test for such problem 
(1). Although the number of robust stability problem is slightly increased, the dimension of the testing 
set is significantly reduced. Furthermore, we considered the uncertainty in inertia matrix M(9d) in 
the above Example, whereas Kamal and Dahleh did not take into account [8]. 



5 Discussions 



The focal point of this paper is to establish a lower dimensional robust stability criterion for robust- 



ness of multivariable control systems. By improving the result given by Kamal and Dahleh, we have 
significantly reduced the high-dimensional stability verification problem to a much lower-dimensional 
problem in parameter space, e.g., for a 2 x 2 interval control system family, our result reduces 8- 
dimensional robust problem to a 2-dimensional one and Dehleh's result reduced this an 8-dimensional 
robust problem to a 4-dimensional one. Furthermore, when dealing with SISO systems, our result 
naturally reduces to the well-known Box Theorem. Of course, from the computational viewpoint, 
checking robust D-stability of interval polynomial matrices is NP-hard. Hence, some sufficient (con- 
servative) conditions were developed based on linear matrix inequality(LMI) techniques. And how to 
reduce the conservativeness of these conditions is still under investigations, on the other hand, from 
the theoretical viewpoint, establishing lower-dimensional necessary and sufficient condition for robust 
D-stability of interval polynomial matrices is important as well. By our result, the 'nucleus' of the 
uncertainty structures can be identified, which is very useful in dealing with robust performance of 
control systems under structured parametric perturbations^' - ^. 

6 Conclusions 

In this paper, robustness of MIMO control systems with interval transfer function matrices is studied. 
For a family of multivariable control systems having interval transfer functions matrices, a lower 
dimensional robust stability criterion has been established. This result improves the previous criterion 
obtained by Kamal and Dahleh, and can be used to establish improved lower dimensional robust 
performance criteria as well. 
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